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Critical Temperature of Ferromagnetic Transition in 
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Yukitoshi Motome and Nobuo Furukawa 1 

Institute of Materials Science, University of Tsukuba, Tsukuba, Ibaraki 305-0006 
1 Department of Physics, Aoyama Gakuin University, Setagaya, Tokyo 157-8572 

(Received ) 

Ferromagnetic transition in three-dimensional double-exchange models is studied by the Monte 
Carlo method. Critical temperature T c is precisely determined by finite-size scaling analysis. 
Strong spin fluctuations in this itinerant system significantly reduce T c from mean-field estimates. 
By choosing appropriate parameters, obtained values of T c quantitatively agree with experiments 
for the ferromagnetic metal regime of (La,Sr)Mn03, which is a typical perovskite manganite 
showing colossal magnetoresistance. This indicates that the double-exchange mechanism alone 
is sufficient to explain T c in this material. Critical exponents are also discussed. 
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ical temperature, finite-size scaling, Monte Carlo method 
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Colossal magnetoresistance near the ferromagnetic 
transition has attracted much attention to the strong in- 
terplay between transport and magnetism in perovskite 
manganese oxides.BEP This feature has been attributed 
to the double-exchange (DE) mechanismEP in the Hamil- 
tonian 

h = - E *« (4<> + h.c.) - j h ]T <n ■ Si, (i) 

ij <7=T4 i 

where Ci a (c\ a ) annihilates (creates) a cr-spin electron at 
site i, and Jh is the Hund's-rule coupling between itiner- 
ant e g electrons with cr, spin and localized ti g moments 
Si. 

Recently this simple DE scenario has been reexam- 
ined from the viewpoint of whether model (Q) quanti- 
tatively describes thermodynamics in a typical material 
Lai_ ;r Sr ;r Mn03 (LSMO), or whether another mechanism 
is necessary. Based on mean-field (MF) arguments, Millis 
et al\y claimed that the DE mechanism alone is insuffi- 
cient to reproduce experiments. In particular, their esti- 
mate for the critical temperature T c is one order of mag- 
nitude higher than experimental results. On the other 
hand, the validity of the DE scenario has_been claimed 
by dynamical mean- field (DMF) studies.EP Magnetism 
as well as transport and optical properties in the LSMO 
compound is well reproduced by model ([!]). 

The DMF method includes dynamical fluctuations 
which were completely neglected in the MF study. DMF 
results suggest the importance of the fluctuations in this 
itinerant system. At the same time, the DMF method is 
not sufficient on its own since it takes into account only 
local fluctuations. In order to settle the controversy of 
whether the DE mechanism alone explains the thermo- 
dynamics of LSMO, it is necessary to make calculations 
which include all fluctuation effects. 

A numerical study is one of the promising candidates 
for including the fluctuation effects appropriately. Many 



surveys based on the Monte Carlo (MC) method have 
been performed. However, their results have not yet con- 
verged mainly due to finite-size effects. Table I summa- 
rizes the values of T c from the MC studies as well as by 
other approximations. MC calculations are performed 
for the system on finite-size clusters. In order to discuss 
critical phenomena, it is necessary to apply finite-size 
scaling analysis. However, such a systematic study has 
not been performed thus far since it has been difficult to 
study large-size clusters in three dimensions (3D) due to 
the rapid increase of the cpu time with the system size. 

In this work we investigate the ferromagnetic transi- 
tion in model (Q) by finite-size_scaling analysis of MC 
results. A new MC algorithmliiP is applied to reduce 
the cpu time markedly, which enables us to study larger 
clusters than ever. Our results elucidate effects of ther- 
mal and spatial fluctuations quantitatively. Our estimate 
for T c agrees with those for LSMO in the ferromagnetic 
metal regime when parameters are chosen appropriately. 
Critical exponents are also examined in comparison with 
experimental results. 

Model (Q) is studied with a classical treatment of the 
localized moments S, as in other surveys listed in Ta- 
ble I. Then the fermion Hamiltonian defined for a given 
configuration of {Si} becomes quadratic in the fermion 
operators and can be diagonalized with the matrix size 
proportional to the system size. Configurations of {Si} 
are stochastically sampled in a MC fashion. The bottle- 
neck of the standard MC techniqueEP is the cost of the 
cpu time required to diagonalize the fermion Hamilto- 
nian in each MC sample. The cost of cpu time is of the 
order of A 3 wiifire A is the system size. We apply a 
new algorithnO which replaces the exact diagortaUza.- 
tion by a moment expansion of the density of states £3l5 
Calculations of the moments can be easily performed on 
parallel computers. The cpu time is reduced to the order 
of (A 2 log A) /Ape where Ape is the number of nodes in 
the parallel calculation. For details of the method, read- 
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ers are referred to rcf. [Tl] . 



Table I. The critical temperature T c of ferromagnetic transition 
in the DE model (Q) estimated by various methods. (MF: mean- 
field approximation, HTE: high-temperature expansion, DMF: 
dynamical mean-field theory, MC: Monte Carlo calculation) The 
results are for x = 0.5 and Jjj = oo except for noted cases. 
Errors in the last digit are shown in parentheses. W is the half 
bandwidth of noninteracting electrons. The third column is T c 
when W is taken as 1 eV. The experimental result (EXP) is also 
denoted for comparison. See text for details. 
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0.0809(3) 
0.0714(l) b 
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0.018(1) 
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0.0208 
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MC 



0.037(3) c 
0.0362(8) d 
0.022(3) b - c 
0.0210(7) b > d 



430(40) c 
422(9) d 
260(40) b ' c 
245(8) b - d 



present 
work 
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a. x = 0.3 and Jh = oo, b. X = 0.5 and Jh = iW 

c. without assumption for the exponent (3 

d. with assumption j3 = 0.365 

e. x = 0.3 



We consider the hopping tij = t only between nearest- 
neighbor sites on a 3D cubic lattice. In the following 
we study the ferromagnetic transition for two param- 
eters, Jh = oo and AW, where W is the half band- 
width of noninteracting electrons (W — 6t in 3D). For 
the case of Jh = oo, components of conduction elec- 
trons with spins antiparallcl to localized moments are 
projected out, and therefore model (|l]) simply takes 
the spinless form of TL = —tJ2ij( c l c j + h- c -)> where 
t = t[cos(9 i /2)cos(9 j /2)+sm(6 i /2)sm(9 J /2)e i ^-^]H> 
The matrix dimension of the Hamiltonian is halved and 
the cpu time is reduced. Moreover, the moment expan- 
sion converges rapidly due to the simple band structure. 
In the following calculations, we perform the moment 
expansion up to the order of M = 20 for Jh = oo and 
M = 40 for Jh = 4H /r _iii order to ensure sufficient preci- 
sion for MC updates. 13 For both values of Jh, we con- 
trol the chemical potential so that the electron density 



) becomes 0.5. Here, the bracket denotes 



the thermal average for the grand canonical ensemble. 
The critical temperature T c in model (Q) is known to 
have a maximum at x ~ 0.5. A closed-shell condition is 
necessary to obtain full polarization in the ground state 
(T = 0). We chose the system size N = N x x N y X N z 
as 4x4x4, 6x6x4, 6x6x8 and 8x8x8 under 
the periodic boundary conditions for the x and y direc- 
tions and the antiperiodic boundary condition for the z 
direction. We have typically run 10, 000 MC samplings 
for measurements after 1,000 MC steps for thermaliza- 
tion. Measurement is divided into five bins to estimate 
the statistical error by the variance among the bins. 
We determine the critical temperature T c by two quan- 



tities. One is the Binder parameterS 1 defined by 

5 =l-(q 4 )/3( 9 2 ) 2 , (2) 

where q is the order parameter which is taken as J^- Si 
in this case. For T > T c (T < T c ), the Binder parame- 
ter decreases (increases) as the system size N increases. 
Thus T c is determined by a crossing point of g for various 
N. The other quantity is the magnetization defined by 

m= [S(k = Q)/N]^ 2 , (3) 

where the spin structure factor S(k) is given by 

S(k) = '£(SiS j )e ik - r ^/N. (4) 

ij 

The thermodynamic limit of m is obtained from the 
system-size extrapolation of S(0) . 

The case of Jh = oo is discussed first. Figure 1 shows 
the Binder parameter for this case. The Binder param- 
eters for different system sizes cross at T ~ 0.03 5 W al- 
though it is difficult to determine T c precisely due to the 
statistical error bars. 

We next study the magnetization. Figure 2(a) displays 
the system-size extrapolation of S(0)/N. The data for 
T < 0.035VF scale well with N~ 2 / 3 , and are extrapolated 
to finite values in the limit of N = oo. This scaling is 
consistent_with fc 2 -dependence of the energy cutoff for 
magnonsEj-* The extrapolated values are summarized as 
the magnetization in Fig. 2(b). The onset of the finite 
magnetization m agrees well with the crossing point of 
the Binder parameter in Fig. 1. The estimate for T c is 
T c ~ (0.037 ± 0.003) 
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Fig. 1. Binder parameter for Jjj 
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To determine the value of T c precisely, we fit the mag- 
netization data by assuming the scaling law of m oc 
(T c — T) 13 . The value of the exponent f3 has not been de- 
termined thus far. At present, our data are insufficient to 
make estimates for both T c and /3 simultaneously within 
the desired accuracy. Instead we assume the exponent in 
3D Heisenberg models (3 = 0.365J13-' and determine the 
value of T c . The data are well fitted, as shown in Fig. 2 
(b). The fit gives T c = (0.0362 ± 0.0008)1^. 

The assumption of (3 in the universality class of 3D 
Heisenberg models is tested by finite-size scaling analysis. 
The scaling hypothesis gives the scaling relation for the 
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spin structure factor (|]) in the form 

S(0) = L 2 -y(L\T-T c n, (5) 

where L — N 1 ^ 3 and / is the scaling function. The crit- 
ical exponents are predicted to-.be 77 = 0.034 and v = 
0.705 in this universality class.O We plot 5(0)/ X 2 -°- 034 
as a function of L\T -T c \ - 705 with T c = 0.0362IF in Fig. 
3. The data appear to follow a universal function (there 
are two branches for T > T c and T < T c ) in this temper- 
ature range. This is consistent with the assumption of 
(3 = 0.365 in the fit in Fig. 2(b). Later we will discuss the 
critical exponents in comparison with the experimental 
results. 
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Fig. 2. (a) System-size extrapolation of the spin structure factor 
and (b) magnetization in the thermodynamic limit for Jjj = 00. 
The curve in the figure is the least-squares-fit to m oc (T — Tc)' 3 . 
See text for details. 
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Fig. 3. Finite-size scaling for Jh = 00. 

We now turn to the case of Jh = AW. As in the 
case of Jh = 00, an estimate for T c is obtained more 
accurately by the magnetization measurement. Figure 4 
shows (a) the system-size extrapolation of S(0)/N and 



(b) the magnetization in the thermodynamic limit. We 
haveT c ~ (0.022±0.003)IF from Fig. 4(b). If we assume 
the scaling m oc (T c — T) 13 with the 3D Heisenberg ex- 
ponent (3 = 0.365 as in Fig. 2(b), the fit for T < 0.02W 
gives T c = (0.0210 ± 0.0007)IF in this case. 
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Fig. 4. (a) System-size extrapolation of the spin structure factor 
and (b) magnetization in the thermodynamic limit for Jjj = 4W. 
The curve in the figure is the least-squares-fit to m oc (T — Tc)' 3 . 
See text for details. 

We compare the results of T c with the other theoreti- 
cal estimates listed in Table I. Our results are one order 
of magnitude lower than the MF estimates. This eluci- 
dates the importance of the strong fluctuations in this 
system. The MF approximation neglects all the fluc- 
tuations, i.e., the thermal and spatial ones. The DMF 
theory gives an intermediate approximation in the sense 
that it includes only the thermal fluctuation and neglects 
the spatial one. Comparison among MF, DMF and our 
MC results clarifies the fluctuation effects quantitatively: 
the thermal fluctuation significantly suppresses T c from 
0.314 7 to ~ 0.0814 7 , and moreover the spatial fluctuation 
reduces T c by more than half from ~ 0.081^ to ~ 0.04M /r 
in the case of Jh = 00. Both fluctuations are important 
and cannot be neglected in a discussion of the phase tran- 
sition. 

The previous MC results listed in Table I deviate con- 
siderably from our results. These deviations arise from 
the finite-size effects and the method used to determine 
the transition temperature. T c has been estimated by 
spin correlation between the faf 
derivative of the magnetizatior 
the data of a finite-size cluster .E 
phase transition in the strict sense, these methods are 
uncontrolled in the absence of system-size extrapolation. 

The values of T c are compared with the experimental 
result in Table I. Here, we take the half bandwidth W 
to be 1 eV, which is typical for 3d transition metal ox- 



thest sites,EP the second 
V and the scaling fit for 
In order to discuss the 
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ides. The bandwidi.h_MZ^= 1 eV is also consistent with 
band calculations Ji§BEJ ) Thus our results, T c ~ 430 K 
for J H = oo and T c ~ 260 K for J H = 4VF, agree well 
with the experimental values of T c in the ferromagnetic 
metal regime, for instance, T c — 369 K at x — 0.3. More 
quantitative agreement will be obtained by tuning the 
value of Jn/W and the doping concentration x. We con- 
clude that once the fluctuation effects are fully taken into 
account, the DE mechanism gives a quantitative descrip- 
tion of T c in LSMO as a major origin. 

Our results strongly support a view of LSMO in the 
ferromagnetic metal regime around x — 0.3 as a canoni- 
cal DE system, which has been pointed out by the DMF 
studies.LF Note that the compound is far from an aii- 
tiferromagnetic/ferromagnetic insulator at x ^ 0.15!!' 
and an >l-type antiferromagnetic metal at i ^ 0.5.EJ 
These instabilities are beyond the scope of model (!]). 
Detailed comparison between experiments involving this 
compound and theories in model (0) provides a good 
starting point to understanding not only the above in- 
stabilities but physical properties in other related mate- 
rials. MC study of other physical properties in 3D is a 
subject for further study. 

In our analysis, the universality class of 3D Heisen- 
berg models is assumed. We also analyze the data by 
assuming the universality class of the mean-field theory 
i.e., (3 = v = 0.5 and rj = 0. In this case, the finite-size 
scaling plot (||) does not show good convergence com- 
pared to the data in Fig. 3. Our data fit better with the 
3D Hciscnberg exponents than do the mean-field ones. 
However, the present data are not sufficiently precise to 
either confirm the universality class of the DE model as 
that of the 3D Heisenberg model, or to identify it as a 
new universality class with similar exponents. 

Concerning the experiments, values of the critical ex- 
ponents for LSMO remain controversial. Ei 
are scattered in the range of 0.3 55 j3 ^ 0.5.E2I 
A recent study claims that the estimate for the criti- 
cal exponent may depend xin the range of the critical 
region assumed for fitting!^ /3 ~ 0.5 is obtained for 
the relatively narrow critical region (T c — T)/T c <; 0.1 
while [3 ~ 0.3 is found for the relatively wide criti- 
cal region (T c — T)/T c <; 0.2. In our calculation, data 
of the magnetization are fitted over the wide range of 
(T c — T)/T c <j 0.4. Considering the experimental situa- 
tion, we cannot exclude the possibility that the results 
for the exponents in the DE model may be modified when 
the system further approaches to the critical point. Fur- 
ther investigation is required to determine the value of 
the exponents with precision. 

To summarize, we have investigated ferromagnetic 
transition in the three-dimensional double-exchange 
model by Monte Carlo calculations. The new algorithm 
based on the moment expansion method has been ap- 
plied to study large clusters. The critical temperature is 
determined precisely by finite-size scaling analysis. Com- 
parison with other theoretical surveys quantitatively elu- 
cidates the importance of thermal and spatial fluctu- 
ations in this itinerant system. The critical tempera- 
ture estimated with appropriate parameters agrees well 
with experiments in the ferromagnetic metal regime of 



(La,Sr)Mn03. The double-exchange mechanism alone is 
sufficient to reproduce the ferromagnetic transition tem- 
perature in this material. The critical exponent in our 
analysis is consistent with experimental trends although 
more detailed study is necessary both theoretically and 
experimentally. 
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